Effective One Body description of tidal effects in inspiralling compact binaries 
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The late part of the gravitational wave signal of binary neutron star inspirals can in principle 
yield crucial information on the nuclear equation of state via its dependence on relativistic tidal 
parameters. In the hope of analytically describing the gravitational wave phasing during the late in- 
spiral (essentially up to contact) we propose an extension of the effective one body (EOB) formalism 
which includes tidal effects. We compare the prediction of this tidal-EOB formalism to recently com- 
puted nonconformally flat quasi-equilibrium circular sequences of binary neutron star systems. Our 
analysis suggests the importance of higher-order (post-Newtonian) corrections to tidal effects, even 
beyond the first post-Newtonian order, and their tendency to significantly increase the "effective 
tidal polarizability" of neutron stars. We compare the EOB predictions to some recently advocated, 
nonresummed, post-Newtonian based ( "Taylor-T4" ) description of the phasing of inspiralling sys- 
tems. This comparison shows the strong sensitivity of the late-inspiral phasing to the choice of the 
analytical model, but raises the hope that a sufficiently accurate numerical-relativity- "calibrated" 
EOB model might give us a reliable handle on the nuclear equation of state. 

PACS numbers: 04.25.Nx, 04.30.-w, 04.40.Dg, 95.30.Sf, 
I. INTRODUCTION 

Some of the prime targets of the currently operating 
network of ground-based detectors of gravitational waves 
(GWs) are the signals emitted by inspiralling and coa- 
lescing compact binaries. Here, "compact binary" refers 
to a binary system made either of two black holes, a 
black hole and a neutron star, or two neutron stars. The 
GW signal emitted by binary black hole (BBH) systems 
has been the subject of intense theoretical studies, based 
cither on analytical methods or on numerical ones. In 
particular, recent progress in the application of the ef- 
fective one body (EOB) approach to BBH systems has 
led to a remarkable agreement between the (analytical) 
EOB predictions and the best current numerical relativ- 
ity results [1, 2] (see also [3]). By contrast, much less 
work has been devoted to the study of the GW signal 
emitted by compact binaries comprising neutron stars: 
either black-holc-neutron-star (BHNS) systems or binary 
neutron-star (BNS) ones. During the inspiral phase (be- 
fore contact), these systems differ from the BBH ones 
by the presence of tidal interactions which affect both 
the dynamics of the inspiral and the emitted waveform. 
During the merger and coalescence phase, the presence 
of neutron stars drastically modifies the GW signal [4- 
6] . The coalescence signal involves (especially in the BNS 
case) a lot of complicated physics and astrophysics, and 
is, probably, not amenable to the type of accurate ana- 
lytical description which worked in the BBH case. Early 
works on this problem have tried to approximately relate 
some qualitative features of the merger GW signal linked, 
e.g., to "tidal disruption", to analytically dcscribablc in- 
puts [7-9]. 

Recently, Flanagan and Hinderer [10-12] have initi- 
ated the program of studying the quantitative influ- 
ence of tidal effects [11, 13, 14] in inspiralling BNS sys- 



tems. However, they only considered the early (lower 
frequency) portion of the GW inspiral signal, mainly be- 
cause they were using a post-Newtonian based descrip- 
tion of the binary dynamics whose validity is restricted to 
low enough frequencies. In particular, one of the results 
of the recent work of Hinderer et al. [12] is to show that 
the accumulated GW phase due to tidal interactions is, 
for most realistic NS models of mass M ~ 1.4M Q smaller 
than the "uncertainty" in the PN-based description of 
GW phasing (see the central panel of their Fig. 4 where 
the thin-dashed and thin-dotted lines are two measures of 
the PN "uncertainty". [These measures are larger than 
the inspiral tidal signal except for the extreme case where 
the radius of the 1.4M Q NS is taken to be > 16 km]. 

By contrast, our aim in this work will be to propose 
a way of describing the binary dynamics (including tidal 
effects) whose validity does not have the limitations of 
PN-based descriptions and therefore is not apriori limited 
to the low frequency part, but extends to significantly 
higher frequencies. This might be crucial to increase the 
dctcctability of the GW signal and thereby have a han- 
dle on the nuclear equation of state (EOS). Indeed, our 
proposal consists in extending the EOB method by in- 
corporating tidal effects in it. Our hope is that such a 
tidally-extended EOB framework will be able to describe 
with sufficient approximation not only the early inspi- 
ral phase, but also the late inspiral up to the moment 
(that we shall consistently determine within our scheme) 
of "contact" . We think that the present EOB description 
of tidal effects is likely to be more accurate than any of 
the possible "post-Newtonian-based" descriptions involv- 
ing supplementary tidal terms (such as [10] or [12]). This 
should be especially true in the BHNS systems which, in 
the limiting case jtins <C m-BH, are known to be well 
described by the EOB approach (and rather badly de- 
scribed by post-Newtonian-based approaches). We will 
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give some evidence of the validity of the EOB description 
of close neutron star systems by comparing our analyt- 
ical predictions to recently calculated quasi-equilibrium 
neutron star (NS) sequences of circular orbits [15] (see 
also [16]). 



II. EFFECTIVE-ACTION DESCRIPTION OF 
TIDAL EFFECTS IN TWO-BODY SYSTEMS 



A. General formalism 

The general relativistic tidal properties of neutron 
stars have been recently studied in Refs. [11-14]. As 
emphasized in [13], there are (at least) three different 
types of tidal responses of a neutron star to an exter- 
nal tidal solicitation, which arc measured by three dif- 
ferent tidal coefficients: (i) a gravito-electric-type coef- 
ficient G/ii = [length] 21+1 measuring the £ th -order mass 
multipolar moment GM ai ^ ae induced in a star by an ex- 
ternal £ th -order gravito-electric tidal field G ai a( ; (ii) 
a gravito-magnetic-type coefficient Gag = [lcngth] 2f+1 
measuring the £ th spin multipole moment GS ai ,,, at in- 
duced in a star by an external £ th -order gravito-magnetic 
tidal field iT ai ... a< ; and (hi) a dimensionless "shape" Love 
number hi measuring the distorsion of the shape of the 
surface of a star by an external ^ th -order gravito-electric 
tidal field. It was found in [13, 14] that all those coeffi- 
cients have a strong sensitivity to the value of the star's 
"compactness" c = GM/cqR (where we denote by Co 
the velocity of light, to be distinguished from the com- 
pactness c). This means, in particular, that the numer- 
ical values of the tidal coefficients of NS's should not 
be evaluated by using Newtonian estimates. Indeed, the 
dimensionless version of m, traditionally denoted as ki 
( "second Love number" ) and defined as 



where 



2k e = (2^-1)!! 



Gju 

R 2l+ 



(1) 



where R denotes the areal radius of the NS, is typi- 
cally three times smaller than its Newtonian counterpart 
(computed from the same equation of state). A similar, 
though less drastic, "quenching" also occurs for the "first 
Love number" hi. In particular, though Newtonian his 
are larger than 1 (and equal to 1 + 2ki, see Eq. (81) of 
[13]), the typical relativistic values of hi are smaller than 
1. This will play a useful role in our analysis below of 
the moment where the tidal distortion of the NS becomes 
too large for continuing to use an analytical approach. 

It was shown in [17, 18] that the motion and radiation 
of two black holes can be described, up to the fifth post- 
Newtonian (5PN) approximation, by an effective action 
of the form 
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point-mass 



E 

A 
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(3) 



is a "skeletonized" description of black holes, as "point 
masses" . To give meaning to the addition of point-mass 
sources to the nonlinear Einstein equations, one needs 
to use a covariant rcgularization method. Refs. [17, 18] 
mainly used Riesz' analytic regularization, but it was al- 
ready mentioned at the time that one could equivalently 
use dimensional rcgularization. The efficiency and con- 
sistency of the latter method was shown by the calcula- 
tions of the dynamics, and radiation, of BBH systems at 
the 3PN level [19-21]. Let us also recall that the lim- 
itation to the 5PN level in Ref. [18] is precisely linked 
to the possible appearance of ambiguities in BBH dy- 
namics appearing at the level where tidal effects start 
entering the picture. Indeed, it is well-known in effective 
field theory that finite-size effects correspond to augment- 
ing the point-mass action 2 by non-minimal (worldlinc) 
couplings involving higher-order derivatives of the field 
[see [22, 23] and Appendix A of Ref. [24]]. More precisely, 
the two tidal effects parametrized by [ii and <ji corre- 
spond to augmenting the leading point-particle effective 
action, (2), (3), by the following nonminimal worldlinc 
couplings 
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Here 1 Gf 



G^.. ae and Rt 



Hi 



ai...a e are the gravito- 
electric and gravito-magnetic "external" tidal gradients 
evaluated along the worldline of the considered star (la- 
belled by A), in the local frames (attached to body A) 
defined in [25]. If needed, they can be reexpressed in 
terms of covariant derivatives of the Ricmann (or Weyl) 
tensor. For instance, using Eq. (3.40) of [25], the leading, 
quadrupolar terms in Eq. (4) read 



AS, 



nonminimal 



E{i^ f ds A £&£*<»> 

A J 

\4 J dSA Baf3 B Aal3 + • • • j- (5) 



where £^ = 



1 B A 

J D at 



L V C 



with 



CjlvaP = \ e wpo cP af) being the dual of the Weyl tensor 
C, and = dz^/ds being the four- velocity along the 
considered worldline. As explained in Appendix A of 
Ref. [24], one can, modulo some suitable "field redefini- 
tions" that do not affect the leading result, indifferently 



S = 



VgR(g) 

UirG 



^point-mass 3 



(2) 



1 We use here the notation of [25], notably for multi- indices L 
a\ , . . . , ag. 



3 



use the Weyl tensor Cap^ or the Ricmann tensor Rap^ 
in evaluating the £ a p and B a p entering Eq. (5). 

The effective-action terms (4), (5) can be used to com- 
pute the various observable effects linked to the relativis- 
tic tidal coefficients [it and at 2 ■ In particular, they imply 
both: (i) additional terms in the dynamics of the con- 
sidered binary system, and (ii) additional terms in the 
gravitational radiation emitted by the considered binary 
system. Both types of additional terms can, in principle, 
be evaluated with any needed relativistic accuracy from 
Eq. (4), i.e. computed either in a "post-Minkowskian" 
(PM) expansion in powers of G/c 2 ,, or (after a further 
re-expansion in powers of 1/co), in a "post-Newtonian" 
(PN) expansion in powers of 1/cg. Let us remark in pass- 
ing that the PM expansion can be conveniently expressed 
in terms of Feynman-likc diagrams, as was explicitly dis- 
cussed (for tensor-scalar gravity) at the 2PN level in [22]. 

Here we shall use the extra terms (4), (5) as a way 
to add to the description of binary black hole systems 
the effects linked to the replacement of one or two of 
the black holes by a neutron star. From this point of 
view, we shall conventionally consider that the tidal co- 
efficients of a black hole vanish: /if 11 = = af u [13, 14]. 
However, as emphasized in [13], more work is needed to 
clarify whether this is exact, i.e. whether the descrip- 
tion of BBH's by an effective action requires or not the 
presence of additional couplings of the type of Eqs. (4), 
(5), as "counter terms" to absorb dimensional regulariza- 
tion poles cx (D — 4) _1 (such poles are indeed linked to 
the possible ambiguities expected to arise at 5PN in the 
point-mass dynamics; see the discussion in Sec. 5 of [18]; 
see also Sec. 7 of [26]) . We leave to future work a clari- 
fication of this subtle issue. 



B. Leading-Order tidal effects in the two-body 
interaction Lagrangian 

Let us first consider the dynamical effects, implied by 
(4) i.e. the tidal contribution to the "Fokker" Lagrangian 
describing the dynamics of two compact bodies after hav- 
ing integrated out the gravitational field, say 



known at the 3PN level. The supplementary term £ tldal 
in Eq. (6) is of the symbolic form (keeping only powers 
of G and 1/cq) 
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tidal 
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Let us start by discussing the leading order contributions 
associated to each tidal coefficient \it or ag. The leading 
term in the contribution linked to fit is simply obtained 
from (4) by inserting the leading-order value of G^, i.e. 
(L = oi . . . at) 



Gi = [d L U ext (x)] A = d£ 



GM 1 



\ZA ~ Z B \ 



(8) 



where B ^ A denotes the companion of body A in the 
considered binary system (A,B — 1,2), and \za — zb\ 
the distance between the two bodies. In addition d£ = 
a , with = d/dz% denotes the differentiation 
with respect to za that appear after taking the limit 
where the field point x tends to za on the worldline of 
body A. Using 
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TAB 



AB 
AB 



(9) 



where n AB 



( Z A z b ) / r AB j r AB = \za~zb\, and where 
the hat denotes a symmetric trace-free (STF) projection, 
and the fact that (see, e.g., Eq. (A25) of [27]) 



AB 11 AB 



l AB " AB 



(21-1)1! 



(10) 



one easily finds that the leading Lagrangian contribution 
proportional to fit reads 



L(q J 



i\ ^point- 



L 



tidal 



(6) 



B\2 



Here, £P° mt - mass (g ; w ) denotes the (time-symmetric) in- 
teraction Lagrangian following from the point-mass ac- 
tion (2) (say after a suitable redefinition of position vari- 
ables to eliminate higher derivatives). It is currently 



2 More precisely, Eq. (4) describes only the effects that are linear 
in tidal deformations (and which preserve parity). If one wished 
to also consider nonlinear tidal effects one should augment the 
quadratic- only terms (5) by higher-order nonminimal worldline 
couplings which arc cubic, quartic, etc. . . in C^ av j3 and its gra- 
dients. The coefficients of such terms would then parametrize 
some nonlinear tidal effects, which have not been considered in 
the linear treatments of Refs. [13, 14]. 



(2i-l)\\ A {GM ) 



2i+2 
AB 



k f G{M B r ^A_ 
r AB 



(11) 



Here we have used (1) to replace Gfxf in terms of the 
dimensionless Love number kf, and of the areal radius 
Ra of the NS. Note that, in a BNS system, one has to 
add two different contributions: L,,a +L,,b. By contrast, 
in a BHNS system one has only L^a if A denotes the NS. 

Let us also evaluate the leading "magnetic-type" con- 
tribution, i.e. the term cx in (6). It is obtained by in- 
serting in (4) the "Newtonian" -level value of the gravito- 
magnetic quadrupolar field H^/ A exterted by body B on 
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body A. This is given by Eq. (6.27a) of [28], namely 
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BA 



TAB 
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where v gA = v B — v A is the relative velocity between B 
and A. A straightforward calculation then yields 



12 a. 



a (GM B ) 2 



' AB 



VAB 
CO 



flAB ■ VAB 

CO 



(13) 

Note that the leading quadrupolar gravito-magnctic 
contribution (13) is smaller than the corresponding 
quadrupolar gravito-electric contribution 



L. 



\,4 {GM 



B\2 



»6 
AB 



by a factor 



l4 



/ V A B \ I n A B ■ VAB 



V c 



C() 



(14) 



(15) 



In terms of the corresponding dimensionless Love num- 
bers J2 (defined in [13]) and &2, the prefactor 8a A /fj, A 
is equal to the dimensionless ratio j 2 /(4fc 2 ). However, 
it was found in [13, 14] that the magnetic Love number 
j2 was much smaller than /c 2 . Typically, for a 7 = 2 
/i-polytrope and a compactness c A ~ 0.15, one has 
32 -0.02, while k 2 ~ 0.1, so that 8cr 2 //i 2 = j 2 /(4fc 2 ) ^ 
—0.05. In other words, the leading gravito-magnetic in- 
teraction (13) is equivalent (say for circular orbits) to a 
1PN fractional correction factor, 1 + a (vab/cq) 2 , modi- 
fying the leading gravito-electric contribution (14), with 
a = 8(T 2 /^2 = .72/(4/22) ~ —0.05. As we shall discuss be- 
low, the 1PN correction to (14), implied by (4), involves 
coefficients a 1PN of order unity. We will therefore, in the 
following, neglect the contribution (13) which represents 
only a small fractional modification to the 1PN correc- 
tion to (14). On the other hand, we shall retain some of 
the higher-degree gravito-electric contributions. Indeed, 
though, for instance, L^a oc l/r AB formally corresponds 
to a 2PN correction to L^a oc l/r AB , its coefficient is 
much larger than that corresponding to an order-unity 
2PN correction to Eq. (14) [see Table I below]. 

Summarizing: the leading-order tidal contributions to 
the two-body interaction Lagrangian arc (from Eq. (11)) 



- tidal 



e>2 



21+2 

"ab 



t>2£+1 

+ k?(M A f%2 

r AB 



(16) 

where kf denotes the £ th dimensionless Love number of 
a NS [11, 13, 14]. Note that the plus sign in Eq. (16) ex- 
presses the fact that the tidal interactions are attractive. 



Structure of subleading (post-Newtonian) 
dynamical tidal effects 



Leaving to future work [29] a detailed computation of 
higher-order relativistic tidal effects, let us indicate their 
general structure. Here, we shall neglect the effects which 
are nonlinear in the worldline couplings fi A of Eq. (4) 
(e.g. effects oc H2H2) f° r two reasons. On the one hand, 
such effects are numerically quite small, even for close 
neutron stars (as we shall check below). On the other 
hand, a fully consistent discussion of such effects requires 
that one considers a more general version of nonminimal 
worldline couplings, involving terms which are cubic (or 
more nonlinear) in the curvature tensor and its covariant 
derivatives. Indeed, it is easily seen that a nonminimal 
coupling which is cubic in G a b ~ £ Q( 3 contributes to the 
dynamics at the same level that a 1PN correction to the 
coupling quadratic in G a bc- 

In the " quadratic- in-curvature" approximation of 
Eq. (4) the part of the tidal interaction which is pro- 
portional to \i A will have the symbolic structure 



B\2 



V ~ n (GM ) 



1 + GM A + GM B 
+ (GM A + GM B ) 



(17) 



where we only indicate the dependence on GM A and 
GM B , leaving out all the coefficients (symbolically re- 
placed by 1), which depend on positions and velocities. 
The presence of an overall factor (GM B ) 2 comes from 
the fact that Gf{z^) in Eq. (4) (which denotes the reg- 
ularized value of some gradient of the curvature tensor 
as the field point x tends to z^(s^) on the worldline of 
M A ) is proportional to GM , so that it is vanishing 
when M B — * 0, i.e. in the limit of a one-body sys- 
tem. [We are considering here a two-body system; in 
the more general case of an A-body system we would 
have G a (za) oc ^ b ^ a GM B ] In a diagrammatic lan- 
guage (see e.g. [22]) the higher-order terms on the right 
hand side (r.h.s.) of Eq. (17) correspond to diagrams 
where, besides having the basic (quadratic in h^ v ) ver- 
tex ha on the A worldline being connected by two gravity 
propagators to two GM B "sources" on the B worldline, 
we also have some further gravity propagators connect- 
ing one of the worldlines either to one of the worldline 
vertices, or to some intermediate "field" vertex. Note 
that the information about the 1PN corrections to both 
gravito-electric (/if) and gravito-magnetic (ai) multipo- 
lar interactions (of any degree £) is contained in the work 
of Damour, Soffel and Xu [28, 30, 31]. We shall dis- 
cuss below the effect of the subleading (post-Newtonian) 
terms in (17) on the EOB description of the dynamics of 
tidally interacting binary systems. 
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III. INCORPORATING DYNAMICAL TIDAL 
EFFECTS IN THE EFFECTIVE ONE-BODY 
(EOB) FORMALISM 

A. General proposal 

The EOB formalism [32-34] replaces the two-body in- 
teraction Lagrangian (or Hamiltonian) by a Hamiltonian, 
of a specific form, which depends only on the relative po- 
sition and momentum of the binary system, say {q,p). 
For a non spinning BBH system, it has been shown that 
its dynamics, up to the 3PN level, can be described by the 
following EOB Hamiltonian (in polar coordinates, within 
the plane of the motion) : 



H 



EOB 



(r,Pr.,P V ) = M\/ 1 + M^eff ~ 1) 



where 



Hefi = 




(18) 



(19) 



Here M = Ma + Mb is the total mass, v = 
Ma Mb /{Ma + Mb) 2 is the symmetric mass ratio and 
Z3 = 2^(4 — 3i/). In addition we are using rescaled dimcn- 
sionless (effective) variables, notably r — vab/GM and 
p v = P v j (GMaMb) , and p rt is canonically conjugated 
to a "tortoise" modification of r [36]. 

A remarkable feature of the EOB formalism is that 
the complicated, original 3PN Hamiltonian (which con- 
tains many corrections to the basic Newtonian Hamilto- 
nian ip 2 + l/r) can be replaced by the simple structure 
(18), (19) whose two crucial ingredients are: (i) a "dou- 
ble square-root" structure H^ob ~ yl + \J P 2 + ■ ■ ■ , and 
(ii) the "condensation" of most of the nonlinear relativis- 
tic gravitational interactions in one function of the (EOB) 
radial variable: the basic "radial potential" A(r). In ad- 
dition, the structure of the function A(r) is quite simple. 
At the 3PN level it is simply equal to 



A (r) = l-2u + 2vu 6 + a 4 vu 4 



(20) 



where a 4 = 94/3 - (41/32)tt 2 , and u=l/r = GM/r AB - 
It was recently found [1] that an excellent description of 
the dynamics of BBH systems is obtained by: (i) aug- 
menting the presently computed terms in the PN expan- 
sion (20) by additional 4PN and 5PN terms, and by (ii) 
Pade-resumming the corresponding 5PN "Taylor" expan- 
sion of the A function. In other words, BBH (or "point 
mass") dynamics is well described by a function of the 
form 

A°(r) = P5 [l - 2u + 2vu 3 + a 4 vu A + a 5 vu 5 + a Q vu e ] , 

(21) 

where denotes an (n, m) Pade approximant. It 
was found in Ref. [1] that a good agreement between 
EOB and numerical relativity binary black hole wave- 
forms is obtained in an extended "banana-like" region in 



the (05, a§) plane approximately extending between the 
points (0,5,0,6) = (0,-20) and (05, ag) = (—36,-1-520). In 
this work we shall select the values 05 = —6.37, a 6 = +50 
which lie within this good region. 

Our proposal for incorporating dynamical tidal effects 
in the EOB formalism consists in preserving the simple 
general structure (18), (19) of the EOB Hamiltonian, but 
to modify the BBH radial potential (21) (which corre- 
sponds to the point-mass action (2)) by augmenting it 
by some "tidal contribution" . In other words the pro- 
posal is to use Eqs. (18), (19) with 



A(r) = A°(r) + A tida \r) 



(22) 



B. 



Incorporating leading order (LO) dynamical 
tidal interactions 



Let us show that, at the leading order (LO), one can 
use a tidal contribution of the form 



A tidal 



(r) 



e>2 



(23) 



with some dimensionless coefficient kJ . 

Indeed, if we keep only the Newtonian approximation 
of the full EOB Hamiltonian (18), (19) (using A(r) = 
1 + A(r) with A(r) = -2 GM/(cl tab) + • • ■ being 1PN 
small as l/cg -> 0) one finds (with \i = M A M B /M) 



# E OB ~Mcl + \^p 2 



1 



nA(r 



o 



(24) 



which exhibits the role of \ [i A(r) as being the interac- 
tion energy. Decomposing A(r) = A°(r) + A tldal (r), and 
remembering that there is a sign reversal between the in- 
teraction energy and the interaction Lagrangian, we see 
that the terms (16) can be converted in a contribution to 
the A(r) potential of the form (23), if the coefficients kJ 
take the values 



2k 



M B ( R A c 2 Q 



1 M A \G{Ma + M b ) 
M a ( R B c 2 



21+1 



2k 



= 2 



M B \G{M a + Me 
M b Mf kf 



2i+l 



{M A 



±M B ) 2e+1 c 2 A e+1 
M A M 2e kf 



{Ma + Mb) 2 ^ 1 cl e+1 



(25) 



In the second form, we have introduced the compactness 
parameters of the stars: ca = GMa/{Ra Cq). It is in- 
teresting to note that the dimensionless tidal parameters 
that enter the EOB dynamics are (when Ma ~ Mb) the 
ratios kf / c 2 A e+1 , rather than the Love numbers kf . Let us 
also note that the velocity of light cq formally appears in 
the numerator of kJ . This is related to the fact that, con- 
trary to the coefficients of the successive powers of u that 
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TABLE I: Tidal properties of BNS and BHNS system. The 
NS model is obtained using the piece-wise polytropic represen- 
tation of EOS SLy and has compactness c = 0.17385. Other 
properties of the model can be found in Table II. 



Model 


q 


,- T 


T 
«3 


T 


BNS 


1 


73.0426 


165.2966 


509.6131 


BHNS 


i 


1.4959 


0.5416 


0.2672 


BHNS 


10 


0.0726 


0.0054 


0.0005 



In particular, we note that, as a function of q, oc 
q/(l + q) 5 vanishes both when q — ► and q — > oo and 
reaches a maximum value when q = Mbh/Mns = 1/4. 
Moreover the maximum value of is larger than the 
value of k£ for a corresponding equal-mass BNS system 
by a factor 4 6 /5 5 = 1.311. We suggest that the numerical 
study of such astrophysically irrelevant BHNS systems 
(with Mbh/Mns ~ 1/4) can be quite useful for improv- 
ing our understanding of tidal interactions in strongly- 
interacting (near contact) regimes. 



enter the BBH EOB potential A°(r) which are (roughly 
speaking) pure numbers of order unity, the coefficients 
kJ entering the tidal contribution ^4 tldal ( r ) w iH tend to 
be much larger than unity (and to increase with £). For 
instance, we shall typically find that = 0(100). This 
numerical difference makes it consistent to add to A°(r) 
(which is known for sure only up to u A terms, i.e. the 
3PN level) additional terms oc u e + u 8 + ■ ■ ■ that would 
formally correspond to 5PN + 7PN + • • • contributions 
if their coefficients were "of order unity" (at least in the 
parametric sense). 

Finally, to illustrate the typical numerical values of the 
EOB tidal parameters we give in Table I the values of 
kJ for three paradigmatic systems, one equal-mass BNS 
and two BHNS of mass ratios q = Mbh /M^s = 4 and 
q = 10. The neutron star model is described with a "re- 
alistic" EOS SLy (with a piece- wise polytropic represen- 
tation, see below) and has the following characteristics: 
mass M = 1.35M0, compactness c = 0.17385, radius 
R = 11.466 km. Note that the main dependence on the 
equation of state (EOS) in kJ (say for the equal-mass 
BNS case) comes from kJ oc R 2 A i+1 . Therefore, if one 
were considering a NS of different radius (because of the 
use of a different EOS) with the same mass, k 2 would be 
approximately given by k 2 ~ 73(i?yi/11.466 km) 5 

One sees in Table I that the dimensionless tidal param- 
eter k 2 is a strongly decreasing function of the mass ratio. 
This is analytically understood by looking at Eq. (25). If 
the label B refers to a black hole (so that kf =0), de- 
noting q = Mbh/Mns = Mb/Ma, we have kJ = (kJ) a 
where 



C. Parametrizing higher-order dynamical tidal 
corrections 

Above we discussed the leading order (LO) contri- 
bution of tidal interactions to the EOB "radial poten- 
tial" A(r). We also discussed the structure of sub- 
leading (post-Newtonian) contributions to tidal inter- 
actions, Eq. (17). Comparing the structure (17) to 
the part of the EOB action linear in A tlda[ , which is 
proportional to the product of A tldal by reduced mass 
/i = M A M B /(M A + M B ), we see that the general struc- 
ture of the tidal contributions to the A(r) potential is 



.4 



tidal 



M J 



M A M B 
GM A GM B 



M B ^ A {GM B ) 2 



1+2 



^ GM A 



GM l 



(27) 



where we invoked dimensional analysis to insert appro- 
priate powers of the (EOB) radial separation r. [Con- 
trary to the action (17) which also depends on velocities 
(and higher-derivatives), the EOB radial potential de- 
pends only on the radius r.] 

In other words, if we separate, for each multipolar or- 
der, the ha and /ib contributions to A tldal , 



^tidal _ j^i 



l>2 



i>2 



we can write 



A"' =A& 1 



u + a 



A(t) 3 

a 3 u . 



(28) 



(29) 



(«!) 



c 2 / +1 (1 + q) 



2£+l ' 



(26) 



Here ca denotes as above the compactness of the NS. 
Therefore, as soon as the mass ratio q is significantly 
larger than one, we see that (nJ) A contains a small factor 
q~ u that suppresses the tidal contribution. As a conse- 
quence, GW-observable tidal effects will be strongly sup- 
pressed in realistic BHNS systems. Note, however, that 
it might be quite useful to compare numerical relativity 
simulations of "artificial" BHNS systems of mass ratio 
q ~ 1 to their EOB description to probe the analyti- 
cal understanding of the late inspiral and plunge phase. 



where 



(30) 



is the part of A^q , Eq. (23), which is linear in fif, or 



, i.e. 



ftp — 2 kf 



M B ( R A 4 



M A \G(M a + M b ) 
Similarly, one will have 



2£+l 



A* = A& 



1 . B(£) . B(£) 2 , B(l) 3 

1 + a 1 u + a 2 u + a 2 u + 



(31) 



(32) 
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The coefficient represents the next to leading order 

(NLO) fractional correction to the leading order A^ e (i.e. 
a 1PN fractional correction), while represents the 

next-to-next to leading order (NNLO) correction (i.e. a 
2PN fractional correction), etc. These coefficients are not 
pure numbers, but rather function of the two dimension- 
less mass ratios 



X A 



M a 



Ma + M B 

M b 
M a + M B 



i-x a . 



(33) 
(34) 



The coefficients entering Eq. (32) are obtained from those 
entering (29) by the interchange of Xa and Xb, i.e. 



a>n (Xa, Xb) = Un yl ' '(Xb,Xa)- The symbolic struc- 
ture (27) would naively suggest that af ^' is a linear 

combination of Xa and Xb and that ce^ is a combina- 
tion of X A , XaXb and X B . However, as the reformula- 
tion of (17) in terms of an EOB potential (27) involves a 
"contact transformation" that depends on the symmet- 
ric mass ratio v = XaXb (see Ref. [32]), the mass-ratio 
dependence of ctn^ might be more complicated. Note 
that, by using the identity Xa + Xb = 1, one can, e.g., 
express a„ in terms of Xa only. [Then a^''' will be 

the same function of Xb than an of Xa-} Note also 
that, if one wishes, one can, for each value of £ factorize 
the total LO terms —Kju 2e+2 , and write 



A tidai = J2 - K ju 2e+2 Af da \ 



where 



with 



l>2 



Af d ^ = l + a[ e) u + a^u 2 



rv£ (-%n rig OL n 

Kn ~\- tig 



(35) 



(36) 



(37) 



Using Eqs. (4.27) and (4.29) of [30], or Eq. (3.33) of [31], 
together with effective action techniques, a recent calcu- 
lation [29] gave the following result for the 1PN coeffi- 
cient of multipolar order £ = 2, a^ 2 \ namely 



A{2) 



2 XA - 



(38) 



More work is needed to determine the higher degree 
and/or higher order coefficients an''' (Xa,^b), and 
thereby the coefficients off entering Eq. (37) . Below, we 
shall focus on the equal-mass case where the coefficients 
Qn become pure numbers. 

Here we shall explore three possible proposals for in- 
cluding higher-order PN corrections in tidal effects. The 
first proposal consists in truncating Eq. (36) at 1PN or- 
der in a straightforward "Taylor" way, i.e. to consider a 



PN correcting factor to the EOB radial potential of the 
form 



^tidal = 1 , -(/) u _ 



(39) 



The second proposal consists in considering a PN cor- 
recting factor which has a "Padc-resummcd" structure, 
i.e. 



1 1 idal 



(40) 



Our third proposal consists in considering a PN correct- 
ing factor which would result from having a "shift" be- 
tween the EOB radial coordinate and the radial coordi- 
nate appearing most naturally in a Newtonian-like tidal 
interaction (cx l/r 2e+2 ). 



A 



tidal 



a\ u 



-(2£+2) 



(41) 



We use here a different notation for the 1PN coefficient, 
a\ , as a reminder that, for instance, when £ = 2, the 
paramctrization (41) corresponds to a 1PN coefficient in 
the paramctrization (39) given by 



r,( 2 ) 



6 a{ 



(2) 



(42) 



IV. COMPARING EOB TO NUMERICAL 
RELATIVITY RESULTS ON "WAVELESS" 
CIRCULAR BINARIES 

The aim of this section is to compare stationary quasi- 
circular configurations of neutron star binaries computed, 
on the one hand, in the analytical framework outlined 
above and, on the other hand, in the numerical frame- 
work recently implemented by Uryu et al. [15] (see 
also [16]). The quantity from both frameworks that we 
shall compare is the binding energy as a function of 
the orbital frequency fl. 



Tidally interacting BNS circular configurations 
in the EOB framework 



1. BNS binding energy in the EOB framework 

As an application of the formalism discussed so far, we 
consider in this section binaries in exactly circular orbits, 
in absence of radiative effects (these will be discussed in 
the following section). 

As the EOB formalism is based on a Hamiltonian de- 
scription of the conservative dynamics, the stable circu- 
lar orbits correspond to minima, with respect to r, of the 
radial potential H^^ l (r,p v ) = H EO B{r,p r , = 0,^). 
Minimizing H^qq 1 ^', p v ) is equivalent to minimizing the 
corresponding effective Hamiltonian H c g, or, its square, 
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2. BNS binding energy in the PN framework 



(HT a r(r,P v ) = A(r) + f J 
= A(u)+p 2 B(u). 



(43) 



Here, we have used the short-hand notation u = 1/r = 
GM/R and B(u) = u 2 A{u). Minimizing (43) with re- 
spect to r (or, equivalently, u), for a given (scaled) total 
angular momentum p v = J tot /GMfi, yields the following 
equation 

A'(u)+ P lB'(u)=0, (44) 

where the prime denotes a u-derivative. This leads to 
the following parametric representation of the squared 
angular momentum: 



(u 2 A(u))> 



(circular orbits), (45) 



where we use the letter j to denote the value of p v 
along the sequence of circular orbits. Inserting this u- 
paramctric representation of j 2 in Eq. (19) defines the 
w-parametric representation of the effective Hamiltonian 
H e g(u). We can then obtain (at least numerically) H e g 
as a function of x by eliminating u between H c g (u) and 
the corresponding w-parametric representation of the fre- 
quency parameter x = (GMfl/c 3 ) 2 ^ 3 obtained by the 
angular Hamilton equation of motion in the circular case 



1 dH ] 



EOB 



MA(u)j{u)u 2 



where iJ rea i denotes the real EOB Hamiltonian 



flEOB = Mi/l + 2l/ LHeff-1 



(46) 



(47) 



In this situation, the binding energy E b of the system is 
simply given by 




l + 2v(H, 



E b (Q) = H EO b - M 



(48) 

where M denotes, as above, the total mass M = Ma + 
Mb of the system, and where one must eliminate u be- 
tween Eq. (46) and Eq. (48) to express the r.h.s. in terms 
of ft. Note that the function E b (fl) depends also on the 
choice of the following parameters: kJ, and af^\ 

Here we shall focus on the equal-mass case, and con- 
sider the dependence of -Eb(ft) only on (k^, , kJ) and 
restrict the paramctrization of 1PN tidal effects to the 
consideration of a single 1PN tidal parameter 6t\ that is 
taken to be the same for the three values of I that we con- 
sider. In addition, we will incorporate 1PN corrections to 
tidal effects in the three aforementioned functional forms, 
Eq. (39)-(41) and contrast their performances. 



We also want to constrast the performance of the EOB 
approach (which represents a rcsummation of the dy- 
namics of the binary system) with the "standard" non- 
resummed PN-bascd description of the binding energy of 
tidally interacting BNS, as used for instance in Ref. [37]. 
The PN-cxpandcd binding energy is written in the form 



point — mass 



(ft) +E tidal (Q), 



(49) 



where 



E 



point — mass (^) 

19 

V - 



~2 X 



27 



/675 

V"6T 



1 



24 
34445 
576 ~ 



205 



12 



155 
~96 l 



v x 



35 
5184 1 



(50) 



is the 3PN accurate post-Newtonian binding energy of 
two point-masses as function of the orbital frequency pa- 
rameter x = (GMfl/c 3 ) 2/3 [19, 38]. The expression of 
the tidal contribution i? tldal (f2) can be obtained for all 
values of the multipolar index £ by noting the following. 
Any (perturbative) power-law radial contribution to the 
interaction Hamiltonian of the form 

SH(r) = (51) 
is easily shown to contribute a corresponding term 



SE b (n) 



i 



(52) 



where it should be noted that the sign of the tidal con- 
tribution flips between the Hamiltonian and the binding 
energy expressed as a function of the orbital frequency 
(ro denoting the Newtonian value of r corresponding to 
a given circular orbit of frequency Q). As a result, we 
have the leading order contribution to the PN-tidal con- 
tribution 



j^tidal 



(0) 



t>2 



\{2l + 2) 



1 



K Jx 2e+2 . 



(53) 



We shall also explore the effect of correcting E£q by a 
fractional 1PN contribution, i.e. to employ a PN tidal 
contribution of the form 



E tldal {x) = (1 + a'^E^ix) 



(54) 



where the (approximate) link with the previously defined 
a\ is 

11 _ 



-ai. 



(55) 



Here the numerical coefficient 11/9 arises as a conse- 
quence of the factor 2n/3 — 1 in the result above (consid- 
ered for n — 6 and n = 7). 
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B. BNS circular configurations in numerical 
relativity 

1. Numerical framework of Uryu et al. 

In a recent paper, Uryu et al. [15] constructed BNS 
systems in quasi-circular orbits by solving numerically 
the full set of Einstein's equations. The important ad- 
vance of this work with respect to previous analyses is 
the fact that Einstein equations are solved for all met- 
ric components, including the nonconformally flat part 
of the spatial metric. This goes beyond the common 
conformally flat approximation that is usually employed 
for the spatial geometry. The conformally flat approxi- 
mation introduced systematic errors which enter the PN 
expansion already at the 2PN level [see the detailed cal- 
culation in the Appendix B of Ref. [35]]. Consistently 
with this analytical argument, it was found in Ref. [15] 
that the difference between conformally flat and noncon- 
formally flat calculations is so large that it can mask 
the effect of tidal interaction for close systems. See, in 
this respect, the location of the conformally flat (IWM) 
binding energy curves in the two upper panels of Fig. 3 
in Ref. [15]. Below we shall however emphasize that the 
nonconformally flat calculations of [15] still introduce sig- 
nificant systematic errors which enter the PN expansion 
at the 3PN level. 

Since the new nonconformally flat results of Uryu et 
al. represent a definitive improvement with respect to 
previous calculations, it is appealing to see to what ex- 
tent these new results agree with existing analytical de- 
scriptions. We extracted from Ref. [15] the six models 
which present the highest computational accuracy. These 
models were obtained by using EOS labelled 2H,HB,2B, 
SLy, FPS and BGN1H1. These labels refer to piecewise 
polytropic EOS. Note that in the case of SLy, FPS and 
BGN1H1 the corresponding piecewise polytropic EOS 
were proposed in Ref. [39] as approximations to origi- 
nal tabulated EOS. In the case of FPS and SLy, this 
implies that the tidal coefficients kg that we have com- 
puted for this work differ (by ~ 20%) from the ones that 
we had previously computed in Ref. [13] that used the 
original tabulated EOS. For example, in the case of a 
neutron star model described by the SLy EOS and hav- 
ing a compactness c = 0.176 (which corresponds to a 
mass of IAMq), we obtain a dimcnsionlcss Love number 
£_(tab) _ o Q7ggg ( w hich is consistent with the first line of 
Table I of Ref. [12]) if we use the tabulated EOS, while we 
obtain k^ poly ^ — 0.09123 if we use the piece-wise poly- 
tropic EOS. Note that the piece-wise polytropic result 
is 18.5% larger than the tabulated one. This suggests 
that one should refine the piece-wise polytropic approx- 
imation to realistic tabulated EOS by incorporating k 2 
within the set of observables that are fitted. 

Among the six EOS that we retain, three, i.e. 2H, 
HB and 2B, use two polytropic intervals, while the other 
three, i.e. SLy, FPS and BGN1H1, use four polytropic 



intervals. We will thus have one dividing density 3 , de- 
noted by po, for 2H, HB and 2B, and three dividing 
densities, (p 0l pi,p 2 ), for SLy, FPS and BGN1H1. Here, 
Po indicates the dividing density between the lower den- 
sity interval that approximates the subnuclear density 
part of the EOS (the crust) and the supernuclear density 
part. The values of (the base- ten logarithm of) po are 
displayed in the first column of Table I. For all EOS, the 
lower density interval ("crust") is approximated by set- 
ting (r , K ) = (1.35692, 3.59389 x 10 13 ), where K (here 
is in cgs units) gives the pressure p in dyn/cm 2 . The other 
dividing densities (for the four-parameter EOS) are fixed 
as pi = 10 14 ' 7 and p 2 = 10 15 . The corresponding adi- 
abatic indices, {Fi,r2,r3}, taken from [15, 39] are also 
given in Table I. For the implementation of the piece- 
wise polytropic EOS we follow the procedure explained 
in Sec. Ill of [39] and in Sec. IID of [15]. 

For each selected EOS, we computed the sequence of 
equilibrium models with the related Love numbers kg up 
to £ = 4. For the compactnesses corresponding to those 
used in [15] we display in Table II the kgs together with 
the values of mass and radius that we obtained from our 
calculation, to check consistency with the corresponding 
values of Table III of [15]. The small differences (at 
the 10 -3 ) level are probably due to the fact that we use 
the finite-digit value of the dividing density po that they 
published. 

2. Subtracting tidal effects from NR data 

Let us start by noting two facts, that can be checked 
from the analytical expressions above, about the depen- 
dence of the binding energy on the tidal parameters kJ: 
i) this dependence is to a very good approximation linear 
and ii) the numerical effect of the strongly dominates 
over that of the higher degree k^'s. For example, if wc 
take the tidal coefficients listed in Table I (which cor- 
respond to the SLy EOS, which yields a radius ~ 11.5 
km for 1.35Mq, which is in the middle of the realistic 
range of NS radii) we find that the tidal contributions to 
the binding energy would reach, if they were extended 
to the maximum frequency that wc shall explore here, 
namely MO max = 0.060, the following values: the k 2 
contribution to E^/M is ~ 3.6 x 10~ 4 ; the contribu- 
tion is smaller than the k 2 by a factor 0.053, and the kJ 
is smaller than the k 2 one by a factor ~ 3.85 x 10~ 3 . 

These two facts allow us to approximately subtract tidal 
effects from NR data. Indeed, if we assume that the 
binding energy computed with a certain equation state 
(EOS)i is approximately given by 

E b {Q; I) « fco(fi) + (4)j/i2(fi) ( 56 ) 



' 5 Here, following the notation of [39] , we use the letter p to denote 
the rest-mass (baryon) density which was denoted by fi in our 
previous work [13] 
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TABLE II: Properties of NS models considered discussed in the numerical analysis of Ref. [15]. The EOS are represented 
as piece-wise-polytropic functions (on four intervals) as proposed in [39, 40]. For the models considered, the present table is 
compatible with Table III of [15]. From left to right, the columns report: the dividing density between the low-density part 
(the crust) and the higher density part of the EOS; the four adiabatic indices for each polytropic interval, {Fo, Fi, T2, Ts}; the 
compactness c = M/R; the NS mass M and the NS radius R; the Love numbers k 2 , &3 and k,4. 



Model 


log(po) 


To 




r 2 


r 3 


M/R 


M 


R 


k 2 


k 3 




2H 


13.847 


1.35692 


3 


3 


3 


0.13097 


1.3507 


15.229 


0.1342 


0.0407 


0.0168 


HB 


14.151 


1.35692 


3 


3 


3 


0.17181 


1.3507 


11.608 


0.0946 


0.0260 


0.0097 


2B 


14.334 


1.35692 


3 


3 


3 


0.20500 


1.3505 


9.728 


0.0686 


0.0174 


0.0059 


SLy 


14.165 


1.35692 


3.005 


2.988 


2.851 


0.17385 


1.3499 


11.466 


0.0928 


0.0254 


0.0095 


FPS 


14.220 


1.35692 


2.985 


2.863 


2.600 


0.18631 


1.3511 


10.709 


0.0805 


0.0214 


0.0077 


BGN1H1 


14.110 


1.35692 


3.258 


1.472 


2.464 


0.15792 


1.3490 


12.614 


0.1059 


0.0307 


0.0120 



we can use the NR data for two different EOS, labelled 
by (/, J) to compute, separately 



h 2 (Cl) 



E b (I)-E b (J) 



(57) 
(58) 



Most importantly we see that Eq. (57) allows us to 
compute from the binding energies of two BNS sequences 
a third binding energy function, ho, which approximately 
represents the binding energy of non tidally interacting 
neutron stars, i.e. the binding energy curve of two point- 
masses. The result of computing the r.h.s. of Eq. (57) for 
five pairs (/, J) of EOS having sufficiently different 's 
is displayed in the left panel of Fig. 1. Two important 
lessons can be drawn from this figure: i) The subtraction 
procedure defined by Eq. (57) is remarkably able to define 
"tidal-free" energy curves that are essentially on top of 
each other; this confirms that our procedure succeeds 
in subtracting out the EOS-dependence of the binding 
energy curves; ii) However, the resulting "universal" ho 
curve still differs significantly both from the EOB point- 
mass curve (black solid line) and the PN point-mass one 
(black dashed line). This second issue will be addressed 
in the next subsection. 

We shall not display here the result of computing the 
/i2 part of the binding energy curve, Eq. (58), because 
it is more sensitive than ho both to numerical noise (in 
the original NR data) and to the presence of higher-order 
tidal PN contributions. Below, we shall address the is- 
sue of determining the tidal contributions to E b with a 
different approach. 



3. Detecting and subtracting systematic errors in NR data 

Here we address the issue ii) mentioned in the previous 
subsection. Indeed, our subtraction procedure has given 
us access to the "universal", EOS-indcpendent part of 



the energy curve ho. However, we have seen that ho still 
significantly differs from the analytical point-mass mod- 
els. We think that the origin of this discrepancy is the 
presence of remaining "systematic errors" in the current 
nonconformally flat approach to BNS systems. Though 
the nonconformally flat integration scheme of Uryu et 
al. is an improvement over previous work, it is how- 
ever still only an approximation to the exact solution 
describing two BNS interacting in a (conservative) "time- 
symmetric" manner (half- retarded- half-advanced) . Here 
we shall only use the data obtained by Ref. [15] called the 
"waveless" approximation. In their approach, "waveless" 
means setting to zero the time-derivative of the confor- 
mal spatial metric (in a certain gauge): dt'jab = 0. As 
the NR gauge is rather similar to the ADM-TT gauge 
used in the 3PN calculation of the interaction Hamilto- 
nian of a two point-mass system in Refs. [19, 41], we can 
see, by looking at the analytical expression of the 3PN- 
accurate ADM Hamiltonian, that neglecting the terms 
containing 7rJ h T ~ dtjab means neglecting some of the 
terms that contribute at the 3PN level. [The simplest of 
these terms being the "kinetic energy" term proportional 
to J d 3 x(irJ b T ) 2 }. This analytical argument suggests that 
the current NR data miss some 3PN contributions, i.e. 
they miss some terms proportional to x 4 in the bind- 
ing energy curve. We are therefore entitled in assum- 
ing that the discrepancy displayed in the left-panel of 
Fig. 1 between the NR ho and the point-mass analyti- 
cal curves is, to leading order, given by an expression of 
the type AE b (fl) = 8x A with an EOS-independent nu- 
merical coefficient 5 that we expect to be of order unity. 
Indeed, the right panel of Fig. (1) exhibits the fact that, 
by subtracting AE b (fl) = Sx 4 , with 5 = 0.8 (see below) 
from all the individual ho curves, we can reach a good vi- 
sual agreement with both analytical point-mass models. 
[Note that the approximate "best-fit" value of S is mainly 
determined by the discrepancy NR/AR on the lower fre- 
quency part of the panel, say for Mfl < 0.035 where the 
contribution to tidal effects is relatively negligible]. 

The remaining differences in this right panel are com- 
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FIG. 1: Comparison between various 5-corrected /lo's (defined in Eq. (57)) and the EOB (resummed, solid line) and 3PN 
(nonresummed, dashed line) point-mass representations of the binding energy. 



patible with the known level of numerical errors in the 
NR data (see Fig. 4 of Ref. [15]). Indeed [15] has used 
the virial theorem to gauge some of the systematic errors 
in their calculation by comparing two measures of the to- 
tal mass of the system (Komar and ADM) . The resulting 
(absolute value) differences in binding energy, say S v Eb 
are in general at the level 10~ 4 M. We used these dif- 
ferences to estimate formal "error-bars" on the various 
energy curves that we use in this work. More precisely, 
in Eb energy curves we add error bars of one-sided am- 
plitude ±7j6Eb, so that the length of the two-sided error 
bars corresponds to the "virial" error. As Fig. 1 con- 
cerns a quantity, /io, defined as a linear combination of 
NR data (see Eq. (57)), we conservatively estimated error 
bars on the ho curve corresponding to the pair 2B-FPS by 
linearly combining in absolute values the corresponding 
individual errors. We use this error bar to gauge the qual- 



ity of the other ho curves (which do not extend as far in 
the high frequency range). This conservative estimate of 
the total error seems appropriate to the present situation 
where the errors are not random, but rather systematic. 
[Note, however, that these "error-bars" seem to be too 
conservative in the lower frequency part of the panels be- 
cause they exceed the "distance" between the ho curves 
and the point-mass models.] Using these error bars we 
can now roughly estimate a range of acceptable values 
of the NR correcting parameter 5. As illustrated in the 
four panels of Fig. 1, the range 0.4<<5<1.2is such that 
the (5-corrected NR-dcduccd ho curves are within "one 
formal sigma" from both point-mass analytical models. 
We shall use this range below to estimate a correspond- 
ing range of probable values of the 1PN tidal parameter 

-(2) 
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4- Least-square analysis: constraining next-to-leading order 
( 1PN) tidal effects from numerical relativity data 

In this subsection we shall firm up the previous analysis 
and make it more quantitative by using a x 2 procedure. 
For each EOS, labelled by index /, we have 20 NR data 
points, Ref. [15] , E b Tyu (x ni ; /) , where the index nj varies 
from one to twenty. We retain in our analysis six EOS; 
I=(2H, HB, 2B, FPS, SLy, BGN1H1). Let us then define 
the following formal x 2 function, measuring the (squared) 
"distance" between NR and EOB: 

X 2 (ai,S) = ^2 

I,n 

(59) 

Here, x = O 2 ' 3 and the index n runs (for each EOS label 
I) over the sample of numerical data from one to twenty, 
so that x 2 contains 120 terms in all. We are interested 
in studying the dependence of x 2 over the two variables 
(r5, (Si). Here 5 denotes the coefficient of a 3PN subtrac- 
tion to NR data of the type that we discussed in the previ- 
ous subsection (as motivated by the neglect of some 3PN 
terms in the "waveless" approximation). As explained 
above, we shall restrict the variation of 6 to the range 
0.4 < 5 < 1.2. For simplicity, we shall actually sample 
this interval through the three values S = (0.4,0.8,1.2). 
On the other hand, the coefficient ot\ parametrizes pos- 
sible next-to-leading order (NLO) 1PN correction to the 
tidal effects. We will use the three different descriptions 
of NLO tidal effects delineated in Eqs. (39)-(41) above. 

We wish to use the least-square method, i.e., mini- 
mizing the EOB-NR "distance" function x 2 (<5i,<5), to 
constrain the values of (a±,5). However, we find that 
X 2 (cti , 5) remains close (on the scale of the NR error bars) 
to its global minimum in a "valley" which extends over a 
significant region of the (<5i,<5) plane. This means that, 
given the present error level in numerical data, we cannot 
meaningfully and simultaneously select preferred values 
for (ai,6). As a substitute, we shall exhibit the sections 
of the x 2 valley that correspond to the three values of 
6 selected visually above in Fig. 1. In other words, we 
now fix S (to one of its three values) in Eq. (59) and 
consider the dependence of x 2 on ai . The resulting one- 
dimensional plots are exhibited in Fig. 2. 

Each panel of Fig. 2 corresponds to a different mod- 
elization of NLO tidal effects: "Taylor" (upper panel,, 
Eq. (39)), "Pade" (middle panel), Eq. (40) and "radial- 
shift" (lower panel), Eq. (41). In addition, each panel 
contains three curves corresponding to the three above- 
selected values of 5: 5 = 0.4 (dash-dot line, right-most 
curve), 6 = 0.8 (solid-line, middle curve), and 5 = 1.2 
(dashed-line, left-most curve). 

Let us start by focussing on the (solid) curves corre- 
sponding to the "central" value of S, 5 = 0.8. We see 
that the preferred values of ol\ that they select (mini- 
mum of the curves) are a\ 7 for the Taylor model, 
ai w 3.5 for the Pade model and L\\ ss 4.5 for the 
"radial-shift" model. This shows that higher order PN 
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FIG. 2: Sections of the function x 2 ( ( 5i,<5) for three values of 
S. The figure displays the corresponding ranges of allowed 
values of Qi. Note that, for all models, the minima are rather 
shallow. 



terms (differently included in the different models) have 
a significant effect on the determination of ot\. Note 
also that when S = 1.2 all the models tend to favour 
a lower value: ol\ ~ 1. The value of x 2 at ai = and 
S = 1.2 is x 2 (0, 1.2) = 5.665 x 10' 7 . This formally corre- 
sponds to an average (squared) "error level" on the indi- 
vidual NR-EOB energy differences summed in x 2 equal 



(< ry >»; I) 4 Ef^{x n -a u I) 
M n) M 
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to • v /x 2 (0,l-2)/120 = 0.687 x 10~ 4 . This level is com- 
parable to the "virial error" on each individual NR data 
point S v Eb/M ~ 10~ 4 . It is therefore reasonable to use 
this level to select a range of values of ot\. Combining 
this range with the range of values of <5's means that, at 
this stage, the range of values of a± that is compatible 
with the NR data is obtained by taking the level surface 
X 2 (cti,8) = x 2 (0, 1.2) as the admissible bottom of the 
"valley" in the (a\,8) plane. This leads to the following 
admissible ranges: < ol\ < 15.7 for the Taylor model; 
< ai < 4.8 for the Pade model; < an < 7.5 for 
the "radial shift" model. It is clear that at this stage 
the fact that (as we have argued above) the NR data are 
"polluted" by some systematic errors (notably linked to 
unaccounted 3PN effects) prevents us from giving very 
significant constraints on the value of ct\. Note in par- 
ticular that the value ct\ = 5/4 = 1.25 which follows ( in 
the equal-mass case) from Eq. (38) is compatible with the 
present NR data (if we allow 8 = 1.2). In this respect, it 
is interesting to note that if we consider a model of the 
form 



^tidal y 



OL\U + Ct2U 



(60) 



with oi\ = 1.25 and compute the corresponding \ 2 for the 
central value 8 = 0.8, we find that x 2 («2,0.8) reaches a 
minimum around «2 ~ 40. In addition the value of the 
minimum of the x 2 is 3.20 x 10 -7 which is slightly bet- 
ter than the performance of the 1PN Taylor model in 
the upper panel of Fig. 2. This shows again that higher 
PN tidal effects can play an important role and that the 
minimima exhibited (for the central value 8 = 0.8) in 
the three panels of Fig. 2 should be viewed as "effective" 
values of ek\. We note in this respect that a situation 
where higher-PN corrections dominate over the 1PN one 
is not at all exceptional. For instance, the 1PN contribu- 
tion to the EOB radial potential A(r) vanishes, its 2PN 
contribution has a rather small coefficient, 2v, while the 
numerical coefficient of the 3PN contribution va^ is quite 
large and significantly modifies the conclusions that one 
might draw from the first two PN contributions. 

Figs. 3 and 4 illustrate the complementary effects of 
ol\ and 8 at the level of the binding energy 

Fig. 3 focuses on the 2B EOS model and contrasts 
(resummed) EOB (top panel) and (nonresummed) PN 
analytical representations of the binding energy. In both 
cases, the NR binding energy is corrected by the same 
amount, that is we assume that 8 takes its central value 
8 = 0.8. We see on this figure that the effect of the 5- 
correction is comparable to that of the added 1PN tidal 
contribution. Note that the value of the parameter 
needed in the PN expanded case (bottom panel) is sig- 
nificantly larger than the ones needed in the EOB case 
(for all ways of modeling 1PN tidal contributions). 

Fig. 4 illustrates the excellent agreement between the 
EOB predictions (here considered for ai = 1.25) and the 
((5-corrccted, with 8 = 1.2) numerical data for all EOS. 
The fact that the \ 2 minima exhibited in Fig. 2 are at a 
comparable level (for all <5's in the range we considered) 
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FIG. 3: 2B EOS: Explicit comparison between various ana- 
lytical representations of the binary binding energy and (cor- 
rected) numerical relativity data. The correction parameter 
is chosen to be S = 0.8. The upper panel refers to EOB 
(resummed) models. The lower panel to PN (nonresummed) 
models. For EOB NLO effects, we use their Pade representa- 
tion, Eq. (40) with ai = 3.5. For the 3PN NLO model, we use 
a[ = 30. 



indicates that a similarly excellent NR/EOB agreement 
would have been obtained all along an extended valley 
in the (cvi,<5) plane. In view of Fig. 3 the same would 
hold for the NR/PN agreement, at the cost, however, of 
using, on average, significantly larger values of a±. 

Summarizing: the recent numerical data of Uryu et 
al. do exhibit the influence of tidal interactions in close 
BNS systems. However, the presence of systematic er- 
rors in the data (due to an imperfect satisfaction -of the 
helical-Killing- vector condition) partially masks the tidal 
interactions and does not allow for a clean determination 
of the coefficients parametrizing tidal effects (and notably 
their 1PN contributions). 

We recommend that new non-conformally flat simula- 
tions be performed for several values of the radius r$ at 
which the hclical-Killing-vector condition is cut off. By 
studying the dependence of the results on ro, it might 
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Uryu et at: 2H EOS 
EOB + tidal NLO : Taylor 
3PN 




FIG. 4: Global comparison between EOB NLO and NR binding energies. We use the values (on, 5) = (1.25,1.2). The 3PN 
point-mass curve is added to guide the eye. 



be possible to extrapolate the results to infinite value of V. INCORPORATING RADIATIVE TIDAL 

r (as used in analytical calculations), and thereby elim- EFFECTS IN THE EOB FORMALISM 

inate the 3PN-level systematic error <5.t 4 . 

Besides the specific Hamiltonian (18), (19), the other 
key ingredients of the EOB formalism are: (i) a specific, 
"factorized" representation of the multipolar waveforms 
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hg, m , and (ii) a resummed estimate of the radiation re- 
action force !F, which must be added to the conservative 
Hamiltonian dynamics (18), (19). In the most recent, and 
seemingly most accurate, version of the EOB formalism 
the radiation reaction is analytically computed in terms 
of the multipolar waveforms. Therefore, it will be enough 
to estimate here the "tidal correction" to the multipolar 
waveforms hi m . Following the "factorization" philosophy 
of Rcfs. [13, 42] we shall look for tidal-correction factors 
fim al = 1 + sucn tnat EOB waveform would 

read 



/•tidal lO 
Jim n in 



(61) 



Here h\ m is the factorized BBH EOB waveform, in- 
troduced in [42], and augmented by two next-to-quasi- 
circular parameters (01,02) in Ref. [1]. [Note, however, 
that, in view of the smallness of the tidal effects on the 
waveform, /^ al — 1 <C 1, it would be equivalent to 
use (as done for the A(r) potential) an additive ansatz: 
ht m = hl n + h%?.] 

In principle, one can use the effective action (4), (5) to 
compute the tidal contributions to the waveform with any 
required rclativistic accuracy (post-Minkowskian and/or 
post-Newtonian) . 

Here, we shall focus on the leading PN-order tidal cor- 
rection to the leading PN waveform, i.e. the I = 2, m = 2 
partial wave /122 ■ This will provide the leading tidal cor- 
rection to the radiation reaction (which is predominantly 
given by a contribution oc |2f2/i22| 2 )- 

In that case, a shortcut for computing the tidal cor- 
rection /22 dal consists in noting that the quadrupo- 
lar gravito-electric contribution in the action (4) cor- 
responds to adding to the energy-momentum ten- 
sor of point-masses an extra contribution A T,K = 

2 g~ x l 2 <5A n0 nminimai S/Sg^, which describes the tidally 
induced quadrupole moment in each body A. At the lead- 
ing "Newtonian" order this means that the quadrupole 
mass moment M^j of the system will be 
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(62) 



where STF denotes a symmetric trace-free projection, 
and where the second term is the tidally induced 
quadrupole moment. Replacing the Newtonian value (8) 
of Gf, (computed using Eq. (9)) yields 
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(63) 



excitation of the internal radial modes: x n — > 0). In addi- 
tion to the explicit tidal modification oc that appears 
in the first factor of Eq. (63), there is an implicit tidal 
effect coming from the fact that the EOB waveform is 
conventionally expressed in terms of the (instantaneous) 
orbital frequency £1 of the binary system. We must then 
eliminate the relative distance tab in Eq. (63) in favor 
of Q. This is done by using the adiabatic (quasi-circular) 
Kepler law. The latter is modified by tidal forces: 
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(64) 
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Differentiating the leading {£ = 2) tidal Lagrangian 
(16), and keeping only the leading (I = 2) term yields a 
modified Kepler law of the form 



Q 2 r 3 
" 'AB 



GM 



M A r\ B M B r\ B 



(65) 



Using (65) to solve tab in terms of f2, and replacing 
the (tidally-corrected) answer in (63) finally leads to a 
quadrupole moment of the form 



Mt. 



rtidal 
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2 - ij 
P r AB n AB 



(66) 



with a tidal-correction factor 



/•tidal 
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= 1 
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= 1 



S'*(£)'( l+ 'E) 



(67) 



The factor /2 2 dal is the I = 2, m — 2 tidal-correction 
factor which was introduced in Eq. (61). It remains, 
however, to eliminate tab in terms of f2, or, as used 
in the waveform of Ref. [42], in terms of the EOB vari- 
able Vq = ro £1 introduced in [43] : at the leading order 



it is enough to use GM j c\r ab 
yields 



vf,(l + 0(l/c 2 )). This 



ftidal 




where fi = Ma Mb /(Ma + Mb) is the reduced mass of 
the binary system, and where we reduced the first ex- 
pression to the center-of-mass frame. Eq. (63) agrees 
with Eq. (7) of [10] (in the limit where one neglects the 



The result (68) agrees (after squaring it) with Eq. (8c) 
of Ref. [10] (in the limit x„ -> 0). 

Summarizing: we propose to incorporate radiative 
tidal effects in the EOB formalism by inserting in the 
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dominant I = 2, m = 2 waveform, a factor of the form 
/ 2 " dal = 1 



G{M A + M B ) 



M A J 



10 

n 



X (1 



(69) 



where we included a possible 1PN correction to radiative 
tidal effects. One then computes a tidal-corrected radi- 
ation reaction by using this corrected waveform in the 
definition of T given in [42] and [13]. In principle the 
(mass-ratio dependent) coefficient /?i can computed an- 
alytically. It can also be "calibrated" by comparing NR 
data of inspiralling BNS systems to the EOB predictions. 



VI. EOB PREDICTIONS FOR THE MOTION 
AND RADIATION OF INSPIRALLING 
COMPACT BINARIES 

Having defined a specific EOB way of incorporating 
tidal effects in the motion and radiation of inspiralling 
compact binaries (BNS or BHNS) let us study the pre- 
dictions made by the resulting tidally-extended EOB for- 
malism. 



Adiabatic inspiral, "last stable orbit", and 
"contact" 



Let us start by considering the adiabatic approxima- 
tion to the inspiral, i.e. the approximation in which the 
inspiral is described as a sequence of circular orbits. In 
this approximation, a key concept is that of the Last Sta- 
ble (circular) Orbit (LSO). We saw above the equation 
determining, in the EOB formalism, the sequence of cir- 
cular orbits, Eq. (44). For large values of p v , and large 
values of r (i.e. small values of u = 1/r), Eq. (44) has a 
unique solution r = 1/u ~ pz, , corresponding to Newto- 
nian circular orbits. However, when p^ decreases (as it 
does along the sequence of inspiralling orbits driven by 
radiation reaction), the sequence of stable circular orbits 
will terminate at certain values tlso = 1/ u lsOj P^ lso 
where there exists a double root of Eq. (44), i.e. a com- 
mon root of Eq. (44) and 



A"(u)+plB"(u) = 0. 



(70) 



The condition determining the radial location of the Last 
Stable Orbit (LSO) is the vanishing of the determinant 



A' B' 
A" B" 



LSO 



A'(u LSO )B"(u LSO ) 
A"(u hSO ) S'(ulso) = 0. 



(71) 



For instance, in the test-mass limit, and in absence of 
tidal corrections, i.e. for A(u) = 1 — 2u. B(u) = 



u 2 A(u) = u 2 -2 u 3 , Eq. (71) reads -4 (1-6 u L so) = 0, so 
that we recover the classic result rxso = I/^lso = 6 (i.e. 
r LSO = 6 GM) for the LSO around a Schwarzschild black 
hole. On the other hand, when inserting in Eq. (71) the 
complete value of the A function, i.e. the sum (22), where 
A°(r;v) is given by Eq. (21), and A tidal (r) by Eq. (23), 
we see that the LSO predicted by the EOB formalism 
will depend both on the symmetric mass ratio v, and on 
the EOB tidal constants kJ , Eq. (25). More precisely, 
these two types of effects (the ^-dependent ones which 
exist already in BBH systems, and the tidal-dependent 
ones which exist only in BHNS and BNS systems) act 
in opposite directions. Indeed, the ^-dependent contri- 
butions tend to make the radial potential A(r) less at- 
tractive (see Eq. (20)), while the tidal ones make A(r) 
more attractive. As a consequence, ^-effects tend to move 
the radial location of the LSO towards smaller values 
( r LSo(^) < 6GM), while tidal effects tend to move tlso 
towards larger values. To avoid gauge effects, it is conve- 
nient to measure the location of the (adiabatic) LSO in 
terms of the corresponding (real) orbital frequency 
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(72) 



Finally, we conclude that the dimensionless orbital fre- 
quency GMil at the LSO is a function of the dimen- 
sionless parameters v, kJ which tends to increase as 
v increases, and to decrease as kJ increases. We have 
seen above that the tidal coefficients kJ generically take 
rather large numerical values, of order «J = 0(100), 
when I = 2, see Table I. However, they enter the A 
function at a higher order in u than the ^-dependent 
effects. As a consequence, the combination of the influ- 
ences of v and kJ, kJ, . . . leads to orbital LSO frequencies 
which arc sometimes larger, and sometimes smaller than 
the "Schwarzschild value" GMfl Schw = 6~ 3 / 2 = 0.06804. 
This is illustrated in Table III which lists the values of 
twice the orbital frequency (corresponding to the adi- 
abatic gravitational wave frequency u>i m for the domi- 
nant mode I = m = 2) for several compactnesses (0.13, 
0.17, 0.17385, 0.5) and for two paradigmatic systems: an 
equal-mass BNS system and a binary black hole system 
(labelled by its formal compactness c = 0.5). Here we 
took the piece-wise polytropic SLy EOS. Note that one 
NS mass is smaller than the "canonical" 1.35Af© so to 
explore a smaller compactness. If needed, one can con- 
vert the dimensionless freqency 2 GMQ in Hz by using 
GM Q = 4.925490947^ s (= 1.476625038 km ) so that the 
conversion factor between uj = GMu> and / — uj/2tt is 



/ = 



2tt GM 



= 32.3125 uj 



Me 
M 



kHz. 



(73) 



We see that, in a BNS system, the LSO frequency is 
smaller than the "Schwarzschild value" 2 GMf2s c h w = 
l/(3-\/6) = 0.136083 for compactness smaller than about 
0.1704. For such system the radius of the LSO is larger 
than the canonical Schwarzschild 6GM. Note, by compar- 
ing BNS to BBH ones, how tidal effects can significantly 
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change the LSO frequency by more than a factor two! 
The results shown in Table III have been computed us- 
ing the leading order, non-PN-correctcd EOB description 
of tidal effects. The inclusion of next-to-leading order ef- 
fects, notably with ai ~ 6, would double the effect of 
tidal interactions at the LSO and would therefore signif- 
icantly affect the numbers listed in the table. 

In some BNS systems the concept of LSO and LSO 
frequency has only a formal meaning because the two 
NS's enter in contact (slightly) before reaching the LSO. 
This is illustrated in Table III which lists also the value of 
(twice) the orbital frequency at the moment of "contact" , 
i.e. when the EOB radial separation R becomes equal 
to the sum of the two (areal) radii Ra + Rb- [We use 
Rb = 2 GMb when the companion is a BH.] Note that it 
is approximately given by the simple analytical formula 



2GMQ. C 



Xa 
ca 



Xb 
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(74) 



This definition of "contact" relies on the use of the EOB 
radial coordinate. As this coordinate is a smooth defor- 
mation of the usual areal coordinate, we think that it is 
a reasonable definition, and we propose here to use the 
EOB description up to the moment when cither the two 
objects enter in contact, or (if it happens earlier) when 
the orbital frequency f2 reaches a maximum. Note also 
that Table III illustrates the possible effect (for c = 0.13) 
of NLO (1PN) tidal contributions. This effect is very 
significant. The second line of the table indicates that, 
when using a "Taylor" model with a\ = 7.0 (which was 
the minimum of y 2 for the central value of 6), the ar- 
rangement of the LSO and touching radius changes. In 
absence of 1PN correction the contact was reached be- 
fore LSO, while with a\ = 7.0 the contact is reached 
after the LSO, which means that the system undergoes 
a short "plunge phase" before entering in contact. 

In addition to the discussion of the frequency at the 
moment of contact (i.e. when R = Ra + Rb) let us also 
consider the dimensionless parameter measuring the tidal 
deformation of the NS labelled A by its companion B 
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At contact, (R = Ra + Rb), this parameter can be 
expressed in terms of the two compactnesses ca = 
GMa/Ra and c B = GM b /Rb as 
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For a symmetric, equal-mass BNS system, we see that, 



upon contact, e™ 1 



1/ 



It was found 



in [13], and briefly recalled above, that the fractional de- 
formation of the NS A is given by the product tA, 
where the "shape" Love number is of order 0.8 for 
a typical NS compactness. This means that, in a sym- 
metric (or near symmetric) BNS system each NS is only 



deformed by about 10% at the moment of contact. This 
motivates our proposal of using the EOB description up 
to the moment of contact. 

In the case of asymmetric BHNS systems (with A la- 
belling the NS and B the BH) we can reach a similar gen- 
eral conclusion by noticing that the dimensionless func- 
tion R 2 A Rb/ (Ra + Rb) 3 (which depends only on the ra- 
tio Ra/Rb) reaches a maximum value of 2 2 /3 3 =4/27 
when Ra = 2 Rb- As a consequence, we have the general 
inequality 



< 



_4_ Cb_ 

27 c A 



(77) 



In the present case, B denotes a BH (with cb = \) so 
that e ™ ntact < 2/(27 c A ) = 0.074074/c A . Upon multipli- 
cation by h$ ~ 0.8 this yields h£ e ™ ntact < 0.06/ca- As 
NS compactnesses are expected to be larger than about 
0.13, we find that the NS in a BHNS system is expected 
to be always deformed by less than 50% up to the mo- 
ment of "contact" with its BH companion. Actually, the 
reasoning above shows that such large deformations are 
only attained when Ra = 2 Rb, i.e. when the mass ratio 
is equal to 



Mb _ c_b_ Rb_ a ^b_ _ 

M A ~ c A Ra~ 
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For typical NS compactnesses ca ~ 0.15, such a mass 
ratio Mb /Ma ~ 1-67 would correspond to a BH of a 
small mass (Mb ~ 2.3 M©). Larger BH masses will lead 
to smaller deformations of the NS. 

Summarizing, the main conclusions of this subsection 
are that: (i) the EOB formalism predicts that the "quasi 
point mass" description can be applied up to contact, 
without the possibility of a disruption of the NS's in a 
well detached state, and (ii) the divide between the sys- 
tems that undergo a plunge before contact and those that 
don't depend strongly both on the compactness and on 
currently uknown higher PN corrections to tidal effects. 

To end this subsection, let us mention that our re- 
sults are robust under the choice of the EOB parame- 
ters 05 and a,6 entering the BBH radial A°(r) potential, 
Eq. (21). The comparison between the currently most so- 
phisticated version of the EOB formalism and the most 
accurate numerical relativity simulations has constrained 
the couple of parameters (05,05) to lie within a rather 
thin banana- like region in the (0.5,0.6) plane. We have 
checked that the results that we present in this paper 
are quite insensitive to the choice of 05 and 06 within 
this "good" region. The default values that we use in 
the present paper are 05 = —6.37, 06 = +50, which lie 
in the "good" region. To illustrate the insensitivity of 
our results to this choice, let us mention that the value 
of twice the orbital frequency at LSO, 2Mfl^Q (05, 05) 
(for an equal mass BNS system and for c = 0.17), changes 
from the value 0.13605, quoted in Table III, to the new 
value 0.13603 for 05 = —4 and 06 = 24 which lie near 
the upper boundary of the good region of paramaters 
discussed in Rcf. [1]. 
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TABLE III: Adiabatic LSO information for BNS and BBH systems. The NS models are built using the piece-wise polytropic 
SLy EOS. From left to rigth, the columns report: the composition of the binary, the compactness c of the objects, the NS 
mass M, the NS radius R, twice the orbital frequency at the adiabatic (EOB) LSO 2GMQ,^^ , the corresponding LSO radius 
r contact /GM, the "contact" frequency 2GAffl contact and the corresponding radial distance r contact /GM. 



System 


c 


M [Mq] 


R [km] 


2GM«ff ( a 1 b 


rLso/GM 


2GM!T2 contact 


r contact j QM 


BNS LO 
BNS LO 
BNS LO 


0.13 
0.17 
0.17385 


1.0050 
1.3205 
1.35 


11.417 
11.470 
11.466 


0.10208 
0.13605 
0.13902 


7.3991 
6.0111 
5.9163 


0.09590 
0.14060 
0.145061 


7.6923 
5.8824 
5.7521 


BNS^iV.o 


0.13 
0.17385 


1.0050 
1.35 


11.417 
11.466 


0.09056 
0.12185 


8.1698 
6.5120 


0.09834 
0.148750 


7.6923 
5.7521 


BBH„ =1/ / 4 
BBH„ = o 


0.5 
0.5 






0.19285 
0.13608 


4.6186 
6.0000 







B. Phasing and waveform from the non-adiabatic 
inspiral of tidally interacting compact binaries 

Let us now consider the motion and radiation of tidally 
interacting binaries predicted by the full EOB formal- 
ism, i.e. beyond the adiabatic approximation. This is 
obtained by integrating the EOB equations of motion 



dr 
~dt 

dPr, 



a(r) 



dm 



EOB 



dp 



R , 



-a(r) 



dH] 



EOB 



01 



dip _ dH EOB 
dt dp v 



dp v 
~dt 



T 

J to 



(79) 



where a(r) = AD 1/2 , H EO B{r,p r ,,P<p) = ^eob/^, with 
Heob defined by Eq. (18) above, and where the (scaled) 
radiation reaction J- v = J-'tp/fJ. is defined in the way in- 
troduced in [13] improved (see Eq. (3) there), i.e. by 
summing over £ and m the adiabatic multipolar partial 
fluxes corresponding to the newly resummed multipolar 
waves hi m (including the tidal correction (68) in h^-i)- 
In addition, we recall that r = R/GM, t = T/GM, 
p v = P^/GMfx, and that the function A(r) is here 
defined as the sum (22). Concerning the other met- 
ric coefficient D~ l (r) (entering the auxiliary function 
a = {A/ B) 1 / 2 = AD -1 ' 2 ) wc replace it by its standard 
resummation (u = 1/r) 



D -1 (r) = l + 6^u 2 + 2(26-3z/)^u 3 



(80) 



The solution of the ODE's (79) is then inserted in 
the newly resummed (and tidally completed) multipolar 



waves hi m to compute the waveform emitted by the inspi- 
ralling compact binary. Here, we shall focus on the I = 2, 
m = 2 dominant asymptotic waveform \iuifi^ 00 (R h^i). 
Scaling it by G/i = GMv and decomposing it in ampli- 
tude and phase, 



R h 



22 



GM 



= A 22 {t) 



3 -i</>22(*) 



(81) 



we can then consider the dominant "metric" gravita- 
tional wave frequency 0J22(t) = d(j>22{t)/dt. [Note that 
all these quantities are dimensionless. In particular 

LV22 = GMUO^ 5 .] 

Up to now we have discussed an extension of the EOB 
formalism which incorporates tidal effects in both the 
motion and the radiation of compact binaries. However, 
it has been advocated [10, 12, 40] to incorporate tidal ef- 
fects as a modification of one of the non-resummed "post- 
Newtonian" -based ways of describing the dynamics of in- 
spiralling binaries. In particular, the recent Ref. [12] uses 
as baseline a time-domain T4-type incorporation of tidal 
effects. To be precise, let us recall that the phasing of 
the T4 approximant is defined by the following ODEs 



2 s 3 / 2 , 



dt 

dx _ 64 
~d~t ~ T 



■x 5 {a^ lm (x)+a^\x)} (82) 



where aJg ylor is the PN expanded expression describing 
pojnt-mass contributions, and where a tldal is given in the 
equal mass case by [10] 



^tidal 



(x) = 26 K 2 X 



(83) 



Here we shall analyze the (metric) GW phase <p22 as 
a function of the corresponding dimensionless frequency 
UJ22 and study the influence on it of tidal effects. More 
precisely, we give here two different comparisons be- 
tween the EOB predictions and the T4 one. In these 
two comparisons, we keep T4 unchanged and defined by 
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FIG. 5: Accumulated GW phase difference (versus GW fre- 
quency W22) between tidal-EOB (quadrupolar) waveforms and 
a Taylor- T4-based PN waveform with (leading order) tidal 
corrections, Eq. (82). Waveforms have been suitably aligned 
(subtracting a relative time and phase shift) at low frequen- 
cies. The circles on the plot indicate, for each curve, the 
dephasing accumulated up to the "contact" frequencies. 



Eq. (82), with a tidal contribution of the leading or- 
der (LO) type (83). On the other hand, we compare 
this tidal-T4 model to two different tidal-EOB models; 
both models use a tidally modified A function, Eq. (22). 
One model (EOB LO ) uses the LO A tidal , Eq. (23), while 
the other one (EOB NLO ) uses the Taylor NLO A tidal , 
Eq. (39), with ai = 7. Here we consider a BNS equal- 
mass system modelled using the 2H EOS with com- 
pactness c = 0.13097, mass M — 1.35M Q and radius 
R = 15.23 km. 

The quantity which is plotted in Fig. 5 is the differ- 
ence A^f 2 OBT4 ^ 22 ) = 0f 2 OBX (u> 22 ) - 2 r 2 4 (^ 22 ). where 
the label X on EOB takes two values, X=LO for the 
leading order model and X=NLO for the next-to-leading 
order model. To compute this quantity we took into 
account possible shifts in both t (i T4 = i EOB + t) and 
4>((j) T4: = EOB + a). We use here the "two- frequency 
pinching" technique of Ref. [44] to fix suitable values 
of the shifts r and a. We use here two pinching fre- 
quencies which are close to 450 Hz. In other words, the 
phase differences displayed in our figure show the phase 
differences accumulated for frequencies between 450 Hz 
and the contact. Though, the figure does not display 
the phase differences below 450 Hz we have checked that 
they remain much smaller than what they become for 
frequencies higher than 450 Hz. 

In Fig. 5, the solid line (black online) dis- 

plays A0™( W22 ) = </>|P BLO ( W22 ) - $i(u 22 ). 
and the dashed line (red online) A</)f 2 OBT4 (o; 22 ) = 

^ 22 ° B (^22) — 22 4 (w 22 ). The two circles on the curves 
indicate the final moments of "contact" . We added 
two vertical (dashed) lines corresponding to 500 Hz and 
1 kHz. 



The main messages that one can draw from this figure 
are: i) the relative dephasing between EOB and T4 (using 
the same tidal model) grows by more than two radians 
up to contact; ii) the inclusion of higher-order PN tidal 
contributions further increases the relative dephasing by 
nearly two radians more. Note that even if one stops the 
evolution around 1 kHz (which is within the sensitivity of 
some possible configurations of Advanced LIGO) the pre- 
viously discussed accumulated dephasings are still larger 
than one radian. This indicates that the GW phasing of 
the ultimate part of the BNS inspiral is very sensitive to 
tidal effects and also very sensitive to their precise ana- 
lytical modelling, including higher-order PN corrections. 
This makes it urgent to do high-accuracy comparisons 
between accurate NR simulations of BNS inspiral and 
EOB models, so as to accurately "calibrate" the EOB 
description of higher-order PN tidal contributions. 



VII. CONCLUSIONS 

We discussed an extension of the EOB formalism which 
includes tidal effects. The hope is that such a "tidal- 
EOB" formalism will be able to go beyond the present 
PN-based proposals whose validity is limited to the early 
(lower-frequency) portion of the GW inspiral signal emit- 
ted by BNS systems. This formalism allows naturally for 
the presence of higher-order PN corrections to the lead- 
ing (Newtonian) effects. We compared tidal-EOB pre- 
dictions to recently computed numerical relativity data 
of quasi-equilibrium circular BNS sequences [15]. We 
showed how to subtract tidal effects from NR data. Even 
after this subtraction, there remains a systematic dif- 
ference between the "point-mass" NR binding energy 
and its EOB (and PN) analytical correspondant . We 
argue that this difference is due to unaccounted 3PN- 
level effects linked to the imperfect satisfaction of the 
helical Killing vector condition (which should be satis- 
fied for physically waveless solutions). We advocate that 
new nonconformally flat simulations be performed for se- 
quences of helical-Killing-vector cut-off radii so as to al- 
low extrapolation to infinite radius. We also suggested to 
study BHNS circular binaries for mass ratios Mbh/-Mns 
of order unity. 

In absence of such physically waveless NR data, we 
propose to subtract from the current data a term 5x 4 
representing a 3PN correction in the binding energy. We 
could then do a least-square analysis to try to mini- 
mize the (squared) "distance" x 2 between NR data and 
tidal-EOB predictions. Our analysis allowed for 1PN- 
corrections to tidal effects parametrized by ai. We 
found that x 2 remains close to its global minimum in 
a flat valley that extends over a significant region of 
the (ai,S) plane. This means that, given the present 
error level in numerical data, we cannot meaningfully 
and simultaneously select preferred values for 6l\ and 
S. Though this analysis is not fully conclusive, it does 
suggest the need of including higher-order PN correction 
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to tidal effects that significantly increase their dynam- 
ical effect. [In other words, the "effective" value, say 
K 2 ff ( u ) = K J (l + + ct2U 2 + . . . ) , which is relevant 
for the late inspiral is significantly larger, by a factor ~ 2, 
than kJ]- These higher-order PN corrections might come 
not from the 1PN level, but from higher PN levels (see in 
particular the end of Sec. IV, where a 2PN completion of 
a recently computed 1PN correction of order unity was 
shown to be fully compatible with current NR data). 

This emphasizes the need both of higher order ana- 
lytical calculations of tidal effects and of high-accuracy 
numerical relativity simulations of inspiralling BNS sys- 
tems. [We note in this respect that it would be useful 
to refine the piece-wise polytropic approximation to re- 
alistic tabulated EOS (used in this paper) by incorpo- 
rating the relativistic Love numbers, notably k2, within 
the set of observables that are fitted]. We argued that 
such a suitably tidally completed EOB formalism will 
be able to describe the dynamics (and GW emission) 
of inspiralling BNS systems essentially up to the con- 
tact of the two neutron stars. We emphasized that, 
though below the dimcnsionlcss (quadrupolar) GW fre- 
quency GMll>22 ~ 0.04 (which corresponds to a frcqcncy 
of 480 Hz for 1.35A/ Q + 1.35M Q system) the present ana- 
lytical knowledge is possibly sufficient for accurately de- 
scribing the system, the GW phasing becomes uncertain 



by a large amount (~ 4 radians) during the late part of 
the inspiral, because of our current lack of secure knowl- 
edge of higher order PN corrections to tidal effects. This 
makes it urgent to do high-accuracy comparisons between 
accurate NR simulation of BNS inspiral and EOB mod- 
els. When the EOB description of higher-PN tidal effects 
is "calibrated" with sufficient accuracy by using such 
EOB/NR comparisons, we think it will be possible to 
use the EOB formalism to extract from Advanced-LIGO 
data some accurate knowledge of the nuclear EOS (via 
the measurement of the crucial parameter kJ). 
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